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Abstract 

We study transport properties of the finite temperature Sakai-Sugimoto model. The 
model represents a holographic dual to 4 + 1 dimensional supersymmetric SU(N C ) 
gauge theory compactified on a circle with anti-periodic boundary conditions for 
fermions, coupled to Nf left-handed quarks and Nf right-handed quarks localized 
at different points on the compact circle. We analytically compute the speed of 
sound and the sound wave attenuation in the quenched approximation. Since con- 
finement /deconffnement (and the chiral symmetry restoration) phase transitions are 
first order in this model, we do not see any signature of these phase transitions in the 
transport properties. 
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1 Introduction 



Recently Sakai and Sugimoto (SS) [H E] introduced a supergravity model realizing 
holographic dual [Hill] to four- dimensional, large N c QCD with massless flavors. Specif- 
ically, they considered five-dimensional SU(N C ) maximally supersymmetric Yang-Mills 
(SYM) theory compactified on a circle with anti-periodic boundary conditions for the 
fermions 1 , and coupled to Nf left-handed quarks and Nf right-handed quarks localized 
at different points on the compact circle. At weak coupling the model can be repre- 
sented by intersecting D4/D8/D8 brane system in type IIA string theory compactified 
on a circle. Coincident N c D4 branes wrap the compactification circle, while two stacks 
(with Nf branes in each) of D8 and D8 branes are localized at different points on this 
compactification circle. At strong coupling, the wrapped D4 branes are replaced with 
an appropriate near horizon geometry, while the D8 and D8 branes are treated in the 
probe approximation 2 . The probe brane approximation is valid in the low-energy limit, 
and as long as Nf <C N c . The probe approximation of the dual holographic description 
corresponds to a quenched approximation for the fundamental quarks on the gauge 
theory side. 

One of the most interesting aspects of the SS model is that it provides a sim- 
ple holographic realization of the nonabelian chiral symmetry breaking. In IE] the 
authors studied finite temperature confinement /deconfinement and chiral symmetry 
restoration in SS model. Rather interestingly, these two phase transitions are not nec- 
essarily simultaneous: for small enough separation of quarks on the circle, it was found 
jS] that there is a phase of the hot SS gauge theory which is deconfined, but with 
a broken chiral symmetry. As the separation between quarks exceeds certain critical 
value, both the deconfinement and the chiral symmetry restoration occur at the same 
temperature. For all range of parameters, all of these phase transitions are of first 
order. Moreover, each phase, even being thermodynamically unfavorable, i.e., having 
a larger free energy, appears to exist at arbitrary temperature 3 . 

In this paper we study transport properties of the hot SS gauge theory plasma 

^uch boundary conditions completely break the supersymmetry and give masses to adjoint 

fermions of the 5d SYM theory. 

2 In other words, the eight-brane backreaction on the D4 brane bulk geometry is neglected. 

3 This should be contrasted with a model with abelian chiral symmetry breaking where the 
chirally symmetric phase is believed to exist only for sufficiently high temperature H3 1101 EH even 
though the chiral symmetry restoration phase transition is expected to be first order. 
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in the quenched approximation. Since the backreaction of the D8 and D8 branes is 
neglected, effectively, we study the hydrodynamics of near-extremal DA branes wrapped 
on a circle with anti-periodic boundary conditions for the fermions. The latter model 
was discussed in ^2] as the first example of the confining theory constructed within 
gauge theory-string theory correspondence. Since the background geometry satisfies 
condition of [T3J E], the shear viscosity 77 of the SS plasma saturates the universal 
viscosity bound proposed 4 in [THj 

V - = -, (LI) 
s An 

where s is the entropy density. On the other hand, the speed of sound and the sound 
wave attenuation is gauge theory specific. Given the dispersion relation for the sound 

waves 

where v a , ( are the plasma sound speed, and bulk viscosity correspondingly, for the SS 
model (in the deconfined phase) we find 

1 C 4 . . 

V5 V 15 

Notice that transport coefficients (jl.3j) are not dissimilar from the transport properties 
of other examples of strongly coupled near-conformal gauge theory plasma [T?| ITH] 

v 2 _i)<i, i~_ K ( V 2 iy K „ lf (L4) 



3/ V V 3 

In fact, the precise value of k = 2 from (jl.3)l is exactly the same as for the cascading 
gauge theory ^H] ■ Also, there is no signature of the confinement / deconfinement phase 
transition in the sound wave dispersion relation. The latter is not unexpected, given 
that this phase transition is a first order. 

The rest of this paper paper is the derivation of (|1.3|) . In the next section we dis- 
cuss effective five- dimensional action of the near-extremal wrapped DA brane system. 
Interestingly, the resulting five- dimensional supergravity action is very similar to the 
one obtained in [TTJ [TH] . In section 3 we study fluctuations of the corresponding black 
brane geometry dual to a sound wave mode of the quenched SS gauge theory plasma. 
We introduce gauge invariant fluctuations and obtain their equations of motion. These 
equations of motion are valid beyond the hydrodynamic approximation, and for arbi- 
trary temperature. In section 4 we derive and solve fluctuation equations analytically 
4 The universality of the shear viscosity in the supergravity approximation was proven in [1311161 IT^| . 
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in the hydrodynamic limit. Imposing Dirichlet condition on the gauge invariant fluc- 
tuations at the boundary of the background black brane geometry determines JH] 
the dispersion relation for the lowest quasinormal frequency (|1.2jl . Finally, using the 
universality result (jl.lj) we can extract (JO}) . 

It would be very interesting to extend our computation beyond the probe approxi- 
mation. We expect that the speed of sound waves and their attenuation would develop 
in this case dependence on the number of fundamental flavors Nf, and the depen- 
dence on the radius of the compactification circle. Relevant discussion of the related 
supergravity background was presented in |2Uj . 



2 Consistent Kaluza-Klein Reduction to 5-dimensions 



In this section we derive the 5-dimensional effective action from the type IIA 10- 
dimensional supergravity action by reduction on S 1 x S A . 

Consider the type IIA supergravity action in the Einstein frame 5 : 



Sua 



2^10 J Mm 



d w x (-G 



(10)\V 2 



R 



(10) 



1 



V„$V^$-e* /2 |F 4 | 2 



(2.1; 



where $ is the dilaton and F4 is the 4-form field strength, and the following metric 
ansatz: 



ds 2 



G j\jdx dec 



e^g^dxW + e 2f e 8w (dS 1 ) 2 + e~ 2w (dS 4 f 



(2.2) 



where the capital Latin indexes (M, N, . . .) run from to 9 and the Greek indexes 
v, . . .) run from to 4. Moreover, the fields / and w as well as $ and F4 do not 
depend on the coordinates of S 1 x S 4 . The 4-form field strength F4 is given by 



Au s i 



(2.3) 



where A is a constant and oj s t is the 4-sphere volume form. With such an ansatz, we 
obtain: 



( _ G (l 0)) l/2 | F4 j 2 = A 2 e -S(f- W ) { _ g) l/2 e -ff (ff4) l/2 



(2.4) 



3 We use conventions of |2] and keep only relevant fields. 
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where g± is the determinant of the metric of the 4-sphere, and the curvature scalar R (10) 
is given by: 



R 



(10) 



e 3 



40 

BP - 20 g^(d,w)(d u w) - -g^(dj)(d v f) 



+ l 2e -2(/-«) j (2.5) 



where R (5) is the curvature scalar with respect to the 5-dimensional metric g^. From 
(|2.4j) and (|2.5|) . and integrating over S 1 x S* 4 in the action (|2.1|) . the effective 5- 
dimensional action follows: 



5, 



2nV 4 

2^io JM 5 



d 5 x (-g) 1/2 



R {5) ~ ^-(df) 2 - 20(dw) 2 - l-(d^) 2 - V 



where V4 is the volume of the 4-sphere and 

<p = A 2 e^e-^ f+8w -12e-f f+2w 



(2.6) 



(2.7) 



Notice that ()2.6|) is very similar to the five dimensional effective action of the cascading 
gauge theory derived in fl] QHl- From (J2.6j) we obtain the following equations of 
motion: 



3 dV 



(2.1 



1 dV 

Dw — = 

40 dw 



(2.9) 



<9$ 



(2.10) 



R% = y djdj + 20 d,wd u w + l - d^d v <$> + l - 9liV V . 
Consider now the following ansatz for the 5-dimensional metric: 



ds\ 



—c\ dt 2 + (? 2 dx 2 + c| dr 2 , 
the correspondent 10-dimensional line element therefore takes the form: 



(2.11) 



(2.12) 



ds 2 



e -ff [_ c 2 df 2 + c 2 dS 2 + c 2 ^ + e 2f+8 W ^2 + g2(/ _ w) ^2 ^ 



A comparison with the finite temperature bulk geometry of the Sakai-Sugimoto model 
in the deconfined phase which is given by |12[ U] : 
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ds 



ss —9s 

+ 

„4 



R 



D4 



Rd4 

15/8 



l\ 2 



-A(r) dt 2 + dx 2 + (dS 1 ) 
dr 2 



4 x2 



R 



D4 



A(r) 

3/4 



+ r 2 (dS 4 ) 



2vrA^ r 



-UJ s i 



(2.14) 



^ = vr^iVy 3 , A(r) = 1- (y)' 



leads to the following identifications: 



/ 

Cl 

c 3 



In what follows we set g s 
(EHoT) are: 



= 1 lnr-A ln ^ 4 , 

1 . 51 . „ 13 . 

= — In q s H in Rr>d-\ Inr 

4 y 80 80 

_2 1 5 i 

= ^ 3 R 2 m re [A(r)]2, 
= Cl [A(r)]"i , 
-gP R 2 m A [A(r)]~5. 
1 and i?D4 = 1 (in this case A 2 



(2.15) 



|), so that the relations 



w 
f 

Cl 

c 2 
c 3 



— Inr , 
10 

13 l 

— Inr , 
80 

5 1 

re [A(r)]2 , 



2 _i 
r 3 [A(r)] 2. 



(2.16) 



We conclude this section with short comments on the thermodynamics of the black 
brane configuration (|2.14j) . The Hawking temperature T is related to the nonextremal- 
ity parameter r\ as 

I ! 7T 

(2.17) 



1 

T 



3r 



1/2 



From (|2.14|) the entropy density s of the black branes is 



s ex r 



5/2 



oc T 



(2.18) 
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The first law of thermodynamics 

-dP = dF = -s dT , 

where P is the pressure and F is the free energy density, then implies that 

-P = F = — sT e = - sT P = - e, (2.19) 
6 6 5 

where e is the energy density. From the equation of state (j2.19J) we find that 

9 dP 1 , 

In the next section we reproduce (|2.2()jl from the dispersion relation for the pole in 
stress-energy tensor two point correlation function in the sound wave channel, or equiv- 
alently from the dispersion relation for the lowest sound channel quasinormal mode 
in the non-extremal geometry ([2.14)1 . 

3 Fluctuations 

Now we study fluctuations in the background geometry 

/->/ + <*/, 

(3.1) 

w — > w + 5w , 
$ -> $ + 5$ , 

where {g^, f, w, <&} are the black brane background configuration (satisfying ([2.16)1 ). 
and {h^, 5f, 5w, 5$} are the fluctuations. We choose the gauge 

ht r = h x . r = h rr = . (3-2) 

Additionally, we assume that all the fluctuations depend only on (t, x^, r), i.e., we have 
an 0(2) rotational symmetry in the X1X2 plane. 

At a linearized level we find that the following sets of fluctuations decouple from 
each other 

\hx1x2} ) 
{fox\x\ h X2X2 } , 

{htxxi h XlX3 } , (3-3) 
\htx2 1 h X2X3 } , 

\htti h aa = h xlXl -\- h X2X2 , hf x ^, h x ^ X3 , 5w, (5$} . 
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The last set of fluctuations is a holographic dual to the sound waves in quenched SS 
gauge theory plasma which is of interest here. Introduce 



hu 


1 


h tt = e - iuJt+iqx3 


4 


-"it j 


h tz 


=4 


h tz = e -*^«X3 


4 


Htz > 


h 


=c\ 


I ___ -iujt+iqxs, 


c 2 
c 2 


11 aa > 


h zz 


=c\ 


^ _ e ~iujt+iqx 3 






Sf 


=e~ 


iujt+iqxs 






8w 


=e~ 










=e~ 


iuit+iqx-z 






•Hi 


haa hzz j Ha 




H a a ~t~ -^zz 



(3.4) 



where -£^ z , if aa , #zz, ^, p} are functions of a radial coordinate only. Expanding 
at a linearized level Eqs. (J2~%J) - (J2~TT|) with Eq. (J3~T|) and Eq. (|3~4]) we find the following 

coupled system of ODE's 



o =K + H' tt 



In 



2 S~l ' 



c 3 



C l V C 2 / 



2 2 &P «9P 



<9w 



9$ 



(3.5) 



=H"z + #L 



In 



Cic 3 



+ ^uqH a 



(3.6) 



0=H': a + H' a 
4 



In 



C 1 C 2 
C 3 



C l V C 2/ 



+ -cl 
3 3 



07* dp ^ 
— ^ + — + — p 

9/ <9u> 9$ 



(3.7) 



=H'' + ii"' 



In 



c l c 2 
C 3 



+ (H' aa -H' tt ) [lnc 2 ]' 



(3. 



2 2 zap 
+ 3 Cs W 



9p 

dw 



dV 



P 
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=T" + T' 



In 



c 3 J 2 c{ 
>,d 2 V v d 2 V _ <9 2 P 



80 



9/ 2 



<9/<9w " (9/9$ 



r 2 

V 2 C 1 



(3.9) 



o =n" + w' 



In 



3" 1 ' 

CiC 2 



c 3 

9 2 P 



+ \w' [H u - HuV + ' 



40 6 \dwdf 



2 

d 2 P 

dw 2 



llu 2 



n 



d 2 V 
dwd& 



P 



,2^1 



n 



(3.10) 



=p" + p' 



111 



£ic| 

c 3 



2 ^ 

-c 2 ( F + 



9 2 P „ <9 2 P 

+ "77T7T P 



d&dw 



5$2 



(3.11) 



where all derivatives are evaluated on the background geometry. Additionally, there 
are three first order constraints associated with the (partially) fixed diffeomorphism 
invariance 



=u H'< + 



In 



C2 



H u )+q [H' tz + 2 



In 



C2 
Cl 



I y /' F + ^w'Q + $'p 



(3.12) 



=Q H' tt 



In* 

Cl 



if« ] + ^fcu H' tz - g F a 



(/ ( y /' .F + AOw'Vl + $'j9 J , (3.13) 



=[ln Cl c 2 ]'^ - [lnc*]' H' a + % ( u 2 H u + 2uq H tz + q 2 % (H tt - H a 



C 2 



+ cj . Z T + f « + % A- (%fF + 40«>W + *</ <3 ' 14> 



9/ dw <9$ 

We explicitly verified that Eqs. ([3.5j) - ([3.1ip are consistent with constraints ()3.12j) - (j3.14j) 
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Introducing the gauge invariant fluctuations 



Zh =4— H tz + 2 id^ — H aa 1 -— — + H tt 

LO \ UJ Z C 2 C2 / LO C 2 



Z f J~ Tj TT" H aa , 

[lncl" 



7 -o _ W H 
w [lncf]' aa ' 



(3.15) 



n c 

and a new radial coordinate 



Z<$> —p — T- — t- if aa , 
[lnc*]' 



x = — , (3.16) 

c 2 

we find from Eqs. (j3.5J) - (j3.11|) . ()3.12|) - (j3.14|) . decoupled set of equations of motion for 

Z's : 

_ d 2 Z g (3g 2 (2x 2 - 1) + 5oo 2 ) dZ H 
dx 2 x (5 a; 2 — q 2 (3 + 2 x 2 )) dx 

4 f (-c; 2 + g 2 x 2 ) (g 2 (3 + 2 x 2 ) - 5 uo 2 ) - 18 g 2 r A x 2 (1 - x 2 ) 5/3 ) 

+ -- Z H (3-17) 

9 (5c; 2 -g 2 (3 + 2x 2 ))(l-x 2 ) 5/3 x 2 r A 

4 g 2 (- 3 g 2 + 5^ 2 ) 
+ W(5u, 2 -g 2 (3 + 2x 2 )) ( 48 ^ + 9 ^ + 52 ^/) ' 



d 2 Z ldz 4 (25 (l-x 2 )(-^ 2 + g 2 x 2 )+ 243r A x 2 (l 



x 2 ) 2/3 



dx 2 x dx 225 r A x 2 (1 - x 2 ) 8/3 

9 

25 (1 — x 2 ) 

(3.18) 

d 2 Z w 1 dZ* 4 ( 2 5 (1 " ("- 2 + fc 2 ) + 162 r A x 2 (1 - x 2 ) 2/3 ) 



dx 2 x dx 225 r A x 2 (l -x 2 ) 8/3 

+ o (12Z/-Z*) 

25(l-x 2 ) 2 V ' ; 



(3.19) 



d 2 Z$ 1 dZ$ 4 



5 (1 - x 2 ) ( — uj 2 + q 2 x 2 ) + 9 r A x 2 (1 - x 2 ) 2/3 ) 



° dx 2 + x dx 45 r A x 2 (l -x 2 ) 8/3 ^* (3.20) 

" 5 (1 -I 2 ) 2 ( ^"^- 
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Further decoupling occurs if we introduce 



k = 48 Z w + 9 Z$ + 52 Z f . (3.21) 

In this case we find 6 : 

Q _ d 2 Z H (3q 2 (2x 2 - 1) + 5iu 2 ) dZ H 
dx 2 x (5u 2 — q 2 (3 + 2 x 2 )) dx 

4 U-u 2 + gV) (g 2 (3 + 2 x 2 ) - 5 u 2 ) - 18 g 2 r A x 2 (1 - x 2 ) 5/3 ) 

+ - ; Zjj (3.22) 

9 (5cj 2 -g 2 (3 + 2a; 2 ))(l-x 2 ) 5/3 x 2 r A 

4 g 2 (-3g 2 + 5iu 2 ) 

+ 15 cu 2 (5cu 2 -g 2 (3 + 2x 2 )) K ' 



d 2 /« 1 ck 4 (u; 2 - q 2 x 2 ) 
0=— + -— + i ^k. 3.23 

dx 2 a; 9r A x 2 (l-x 2 ) 7 

4 Hydrodynamic limit 

We study now physical fluctuation equations (j3.22j) . ()3.23j) in the hydrodynamics ap- 
proximation, u) — > 0, g — > with - kept constant. Similar to the computations in 
[T7| ITS] , we would need only leading and next-to- leading (in q) solution of ()3.22|) and 
()3.23|) . We find that at the horizon, x — > 0+, Zh oc x ±tu} ^ 2nT \ and similarly for k. 
Incoming boundary conditions on all physical modes implies that 

Z H ( X ) = x' iw z H (x) , k{x) = x~ ivo lC(x) , (4.1) 

where {zh, fC} are regular at the horizon; we further introduced 

'-Sr- (4 ' 2) 

There is a single integration constant for these physical modes, namely, the overall 
scale. Without the loss of generality the latter can be fixed as 



z H {x 



= 1. (4.3) 



6 This set of gauge invariant fluctuations will be sufficient to determine the sound wave dispersion 
relation. 
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In this case, the pole dispersion relation is simply determined as [T9"] 



z H (x 



(4.4) 



The other boundary condition (besides regularity at the horizon and ([4.4)1 ) is [T9"] 



/C(x) 



(4.5) 



X->1_ 



Let's introduce 



2# = z#,o + i q Zff,i , /C = /C + i q /Ci 



(4.6) 



where the index refers to either the leading, oc q°, or to the next-to- leading, oc q 1 , 
order in the hydrodynamic approximation. Additionally, as we are interested in the 
hydrodynamic pole dispersion relation in the stress-energy correlation functions, we 
find it convenient to parameterize 



ro = v s q — i q 2 T 



(4.7) 



where the speed of sound v s and the sound wave attenuation T are to be determined 
from the pole dispersion relation (|4.4jl 







0. 



(4.8) 



Using parameterizations ([4.6)1 . ()4.7[) . we obtain from ([3.22)1 and ([3.23)1 the following 
ODE's 



=x K% + K.' 



(4.9) 



„ (6x 2 -3 + 5^ 2 ) , 8 4(-3 + 5^ 2 )/C 

ZH '° x(2x 2 -5v s 2 + 3) Zh '° 2x 2 -5v. 2 + 3 Zh '° 15 (2x 2 - hv 2 + 3) v 2 ' 

(4.10) 

describing leading (oc q°), and 

U=xK!{ + K! 1 -2v a K! Qi (4.11) 
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„ (6x 2 -3 + 5^ 2 ) , 8 

x (2x 2 -5tv» + 3) 2x2 -5^ + 3 ^ 

2w s (40 x 2 r + 20 x 2 v s 2 - 25 + 30 v s 2 - 4 x 4 - 12 x 2 - 9) , 
+ x(2x 2 -5v s 2 + 3) 2 

8w s (-2 x 2 + 5 v s 2 - 3 + 10 T) 8 T (6 x 2 + 9 + 25 t» s 4 - 30 w s 2 ) /C 

72 ZH,0 + TT 



(2x 2 -5^ 2 + 3)' ' 15 (2x 2 -5^ 2 + 3) 2 ^ 3 

(-3 + 5^ 2 )/d 



15 (2x 2 -5t; s 2 + 3)w s 2 ' 

(4.12) 

describing next-to- leading (oc q 1 ) order in the hydrodynamic approximation. Solving 
(14. 9|) and (|4.11|) subject to regularity at the horizon and the boundary condition (|4.5jl 
we find 

/C = 0, /Ci = 0. (4.13) 

Given (I4.13jl . solution to ()4.9|) subject to regularity at the horizon and the boundary 
condition ()4.3|) is 



5w s 2 + 2x 2 -3 
5 v s 2 - 3 

which from (|4.8jl determines (in agreement with ()2.2()j) 



••=^- (415) 

With ()4.13j) - (j4.15j) . solution to (|4.12jl subject to regularity at the horizon is 

z H ,i = C (l-x 2 ) + 4= (5r-2) , (4.16) 
V5 

where C is an arbitrary integration constant. From ()4.8j) we conclude 7 

r = |. (4.17) 

Finally, comparing ()1.2|) and ()4.7j) . and using fjl.ip . we obtain from ()4.17|) result quoted 
in (USD. 

7 Boundary condition (|4.3() fixes C = 0. 
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